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Abstract — There are many technique to approximate the 
dimension of fractal sets. A famous technique to approximate 
the fractal dimension is correlation dimension. In this paper, we 
propose a new approach to compute the correlation dimension of 
fractals generated by the Escape Time Algorithm (ETA) and 
computing the correlation function by Grassberger-Procaccia 
method, was implemented using the Matlab program. A log -log 
graph is plotted by matlab program, the result will be an 
approximate polynomial of degree one (straight line) that’s fits 
the data in a least square method. The correlation dimension of 
fractal object is the slope of this straight line. 

Index Terms —ETA, geometry, GPA, IFS. 


I. INTRODUCTION 

The fractal geometry is a kind of geometry different from 
Euclidean geometry. The main difference is in the notation of 
dimension. The Euclidean geometry deals with sets which 
exists in integer dimension while the fractal geometry deals 
with sets in integer dimension. Dimension is the number of 
degrees of freedom available for a movement in a space. In 
common usage, the dimension of an object is the 
measurements that define its shape and size. The dimension, 
in physic, an expression of the character of a derived quantity 
in relation to fundamental quantities, without regard to its 
numerical value, while in mathematics, the dimension 
measurements are used to quantify the space filling properties 
of set [5]. Fractal geometry is a useful way to describe and 
characterize complex shapes and surfaces. The idea of fractal 
geometry was originally derived by Mandelbrot in 1967 to 
describe self-similar geometric figures such as the Von Koch 
curve. Fractal dimension D is a key quantity in fractal 
geometry. The D value can be a non-integer and can be used 
as an indicator of the complexity of the curves and surfaces. 
For a self-similar figure, it can be decomposed into N small 
parts, where each part is a reduced copy of the original figure 
by ratio r. The D of a self-similar figure, it can be defied as 
D = — log(A r ) / log(V'). For curves and images that are not 
self-similar, there exist numerous empirical methods to 
compute D . Results from different estimators were found to 
differ from each other. This is partly due to the elusive 
definition of D , i.e., the Hausdorff-Besicovitch dimension 
[4]. There are many kinds of dimensions. The correlation 
dimension is one of the fractal dimension measurement 
because it permits non-integer values. In (1983) Grassbergr 
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and Procaccia [7, 8] published an algorithm (GPA) for 
estimating the correlation dimension that has become very 
popular and has been widely used since. 

The correlation dimension O €or can be calculated in real 
time as the fractal generated by Escape Time Algorithm by 
using the distances between every pair of points in the 
attractor set of N number of points, it is so-called correlation 
function C (e) . It is defined as the probability that two 
arbitrary points on the fractal are closer together than the 
sides of size £ of the cells which cover the fractal. 

The paper is organized as follows: In section 2, some 
background material is included to assist readers less familiar 
with detailed to consider. Section 3 introduces the concepts 
needed to understand the correlation function for calculating 
the correlation dimension D c ^ r , In section 4, we presents the 
Escape Time Algorithm. In section 5, the proposed method to 
compute the correlation dimension of fractals generated by 
Escape Time Algorithm. Finally, the conclusions are drawn in 
section 6. 

II. Theoretical B ackground 

This section presents an overview of the major concepts and 
results of fractals dimension and the Iterated Function System 
(IFS). A more detailed review of the topics in this section are 
as in [1,3,6]. Fractal is defined as the attractor of mutually 
recursive function called Iterated Function System (IFS). 
Such systems consist of sets of equations, which represent a 
rotation, translation and, scaling. Given a metric space (X d ), 
we can construct a new metric space (Ti'OO.. fi ) , where K(X) 
is the collection of nonempty compact set of X , and h is the 
Hausdorff metric on defined by: 
h(A,B ) = max{d(A B), dfB, A)}. 

Definition 1 [9] A transformation / : X X on a metric space 
( X , d " i is said to be a contraction mapping if for some 
O < ?" < 1 difbdjiy')) < rd (x, y), for .all x, y E X, 
Theorem 1 (contraction mapping theorem)^ 3]. Fet /: X -* X 
be a contraction mapping on a complete metric space (A r , d). 
Then / admits a unique fixed point x L in X and moreover for 
any point xeX , the sequence j/ [n] GO : n= 0,1,...] 
converges to x“, where denotes the n-fold composition of 
/. 

Definition 2 [6]. An IFS consists of a 

complete metric space (X d j and a finite set of contractive 
transformation : X ^ X with contractivity factor r n , for 
n = 1 , 2, , N . The contractivity factor for the IFS is the 
maximum r among }. The attractor of the IFS is the 

unique fixed point in Jf(X) of the transformation 

IV: 00 00 . defined by A = IV00 = for 

any B £ J4" 00 . 

Definition 3 [6]. A transformation/:^ -> A r on a metric 
space (X c-j is called dynamical system. It is denoted by 
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{X /} . The orbit of a point x t X is the sequence 

Theorem 2 [ 6 ]. Let }be an IFS with 

attractor A . The IFS is totally disconnected if and only if 
co; (A) n iiij-j- (A) = 0 , for all i,j E [l , 2, , A : } with i == j. 

Definition 4 [ 6 ]. Let { a; = 1, 2, ... , A 7 } be totally 
disconnected IFS with_attractor A , The transformation 
f:A A, defined by /(ft) = , for ft E tt n (.4), is 

called the associated shift transformation on ,4 . They 
dynamical system £4;/} is called the shift dynamical system 
associated with IFS. 


III. Correlation dimension 


There are many ways to define the fractal dimension, but 
one of the numerically simplest and most widly used is the 
correlation dimension of Grassberger and Porcaccia (1983) 
[7,8]. The real interest of the correlation dimension is 
determining the dimensions of the fractal attractor. Ther are 
others methods of measuring dimension such that hausdorff 
dimension, the box- counting dimension but the correlation 
dimension has the advantage of being straight forwardly and 
quickly claculated, when only a small number of points is a 
vailable, and is overwhelmingly concord with other 
calculations dimension. 

Definition 5: Let {* J'rli be a sequence of points in R n . The 
correlation function, C(s) is defined by: 

C (V) = lim JV _ >3B 77 X {number of pairs (x^xjof 
points with d(x^xj~) < er} 

C fe) = liniy^ V £f; 0(E - \\x { - XjH3 

where S(rj is the Heaviside step function (unit step function) 
which has a value of either 0 or 1 and may be defined as: 




:i, o < eo-|l* E -*rll) 
0, O>0( £ - 11*1-^11) 


which acts as a counter of the number of pairs of points with 
Euclidean separation - distance between two points on the 

fractal, x> and x.. The multiplier 77 is included to normalize 


the count by the number of pairs of points on the fractal. . 

Grassberger and Procaccica [7, 8 ] and Barnsley [ 6 ] 
established that for small values of the separation distance z, 
the correlation function C fe) has been found to follow a 
power law such that 
C (f)^Cc~ D - T 

where C (V) is the number of cells with the edge size 

£ necessary to cover all points of the fractal object and C is a 

positive constant. The ' w ' symbol in this expression is used to 

indicate that this is not an exact equality but is a scaling 

relation that expected to be valid for sufficiently large N and 

small £ . After taking logarithms of each side of the scaling 

relation and rearranging terms, so we define 

log C (f) 

D tcr - li"n-— 

log(£Tj 

The correlation dimension D cor estimated using the least 
squares linear regression of log C (/) and versus log(ff), then 
the slope of linear model represent D Cl0r . 


IV. The Escape Time Algorithm [6] 

In this section, we study fractals generated by means of 
Escape time algorithms (ETA). Such fractals are computed by 
repeatedly applying a transformation to a given initial point in 


the plane, the algorithm is based on the number of iterations 
necessary to determine whether the orbit sequence tends to 
infinity or not that is, an orbit diverges when its points grow, 
further apart without bounds. Escape time algorithm is 
numerical computer graphical experiment to compare the 
number of iterations for the orbits of different points to escape 
from ball of large radius, centered at the origin. A fractal can 
then defined as the set whose orbit does not diverge. 

The Escape Time Algorithm applied to any dynamical 
system of the form {M z ;/},{£;/}, or {£;/}. Consider a 
window TV and a region^, to which orbits of points in TV 
might escape. The result will be a "picture" of IV, where in 
the pixel corresponding to the point z is colored according to 
the smallest value of the positive integer n , such that 
/ [n] (z) E 31. Let (ft, b j and (c, ftj be the coordinates of the 
lower left corner and the upper right corner respectively of a 
closed , filled rectangle IV c R L . 

Let an array of points in IV defined by: 

X 04r = (o+p b+q ^^3 , for p,q= 0,1 . M 

where Af is a positive integer. 

Define: 3£ = [(x, y) E R L : x* + y 2 >r}, where r is a 
positive number 

compute a finite set of points: 

1:1 belonging to the orbit 

of x u E IV, for each V,q = 0. 1, ...,3: . The total number 
of points computed on an orbit is at most A : (where A : is a 
positive integer). If the set of computed points of the orbit x p ^ 
contain no a point in 31 when n — N , then the pixel 
corresponding to x 3 ^ will be black color, and the 

computation passes to the next value of (p, g j. Otherwise the 
pixel corresponding to x v is assigned a color indexed by the 

first integer n, such that f^(x pq ) E 31 , and then the 

computation passes to the next value of [p , g j. 

So, the fractal generated by calculating the orbit for each point 
in the plane (that is pixel on screen) and checking whether it 
diverges. In our work, the black and white image of the fractal 
generating by coloring a pixel white if it is the orbit diverges 

i.e., / [f!] (x pifl ) E X, for some n < A r , and black if it is orbit 

does not i.e., (_t l1 ) e 31 for all n < N. So, we have a 

closed subset F of IV constructed by the Escape Time 
Algorithm, defined by: 

F — [.Vp , E TV, such that ; (x p £ 31 for all n < A r } 

= E TV, such that . is black point}. 

The set F is called the Escape Time Fractal. 

The Escape Time Algorithm [ 6 ]: 

1. Given TV Q Q£ z , s.t 

TV = (U, >0 : ft < x < £ 7 , b < y < d) , define the 
array of point in TV by , x pz = (a p(c - a)/M , 
faEq(ft - 

p,q = 1,2 , ...., M, for any positive integer M. 

2. Let C be a circle centered at the origin, the set 31 is 

defined such that, 

J2 = {(x, y) E M : : x 2 + y 2 > r} , where r is 

sufficiently large number. 

3. Let / be a complex function with the orbit of point 

(A 1 where e W. 

4. Repeat, V x p ..~ E TV. 
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IF ■|[/ ,[i1 ' C-X'ti ,)}- =c E 31, then x p .. is colored with the 
color indexed by n. 

Else it is colored Black. 

End IF. 

5. Change all that is not black color to the white color, 
such that: 

if f^(xp -) E 31 for somen < N 
p [F hack if/ [nl (s p ^) £ if for all n < N 

6. Then the set of escape time point -4 in W is defined 
as follows, 

F = {x M E IV : / M (x p ^) e 31 for all n < N} 

F = {.x L1 E W , such that x». : is black point}. 

Output: The set F is called the fractal constructed by Escape 
Time Algorithm. 

To apply the Escape Time Algorithm for IFS we find a 
relationship between the dynamical system {R 2 ,/} and IFS 
[ER 2 ; which is the shift dynamical system 

{F f f] associated with the IFS. 

So, we must follow these steps: 

1. Find a dynamical system { M 2 }/} which is an 

extension of a shift dynamical system associated 
with the IFS, and which tends to transform points 
off the fractal of the IFS to new points that are 
further away from the fractal (this is always 
possible if the IFS is totally disconnected ). 

2. Apply the Escape Time Algorithm, with 3i and IV 

chosen properly. 

V. The Correlation Dimension of Fractal Generated 

by Escape Time Algorithm 

In this section, we propose on approach algorithm to 
compute correlation dimension of fractals generated by 
Escape Time Algorithm by finding shift dynamical system 
associated with IFS and the correlation function has been 
investigated by Grassbger-Procaccia Algorithm[7,8]. 

A log-log graph of the correlation function C (eQ versus the 
distances between every pair of points ^ in the fractals 
constructed by Escape Time Algorithm is an approximate of 
the correlation dimension, i.e., the correlation dimension is 
deduced from the slope of the straight line scaling grid in a 
plot of log C(if') versus lag(f) . We presented and 
implemented in the Matlab program listed in Appendix for 
computing the correlation dimension of fractals generated by 
Escape Time Algorithm. 

The Proposed Algorithm (1): 

INPUT: READ m , i- { m = number of regions, v = 
value } 

FOR f = 1 TO m 

READ 1 (x, y) = ($; (xj , h; (y) ') 

ENDFOR { i } 

READ a , b , c , d , r, M , numits 
PROCESS: FOR P = 1 TO M 

FOR q = 1 TO M 

x = a E (c — □) *p SM 
y = b E (□’ - b) * q /M 
FOR n=l TO nmits 
IF (x {< >,<,>} v) and (y{< >, 
<, >} v) THEN 

* = Si 00 , y = h (y) 

ENDIF 


IF x* x + y*y>r THEN 
Graph(jJ, g), Color (n), n = numits 

ENDFOR { n} 

ENDFOR{ q } 

ENDFOR {p } 

OUTPUT: GRAPH of m Regions 

END. 


{.Calculate correlation dimension (CorD) 

.} 

Stepl: INPUT Nt, N {Nt=No. of transients, N=NO. of 
points} 

Step2: Generate random x ;V[ , i = I : N £ 

Step3: d = xj) 2 F (y £ - yj) 2 , 

= mirj(cj , = m ax{d) 


log 1 


IctetCsI 


El 


Step4: r = 53SBB!Sl (fl = 

log (23 

Step5: = 2(.mx z ) C —i - lXt = 1, 

Step6: A r ti = « 




iVj, = d < t d > 0] ,5 = y Aj. , c£{ = ■connate \ce;, 
] , i = 1: n E 

Step7: Op = 3.. /Ci = E 1 , k 2 — 


xa-i = 


leg :;-e L’jj 


y*i = 


log(2) 


lorgps] 

Xp = xd L t y p = yd I , i = k 1 :: k 2 

Step8: cp E = Poly it „ y D f l) , C orD = 

OUTPUT: PRINT (CorD) 


A -' LI 


In the following, this algorithm have been applied to the 
examples of fractals generated by Escape Time Algorithm. 
That is, we calculated the correlation dimension of fractal 
generated by Escape Time Algorithm. 

Example 1: Consider the IFS co 2 , co 2 }, where 

&i(x,y) = + ^y),at a (x,y} = y) 

The attractor of the IFS is the fractal Sierpinski triangle S. 
The relationship between the dynamical system (R . /} and 
fractals IFS co^. co 2 , ii> g } is a shift dynamical system 
[5,/} associated with IFS. Calculate the inverse 
transformations for the above system, which are as follows: 
Define / » K 2 , by 

f(x,y) = 

(£rjf 1 (.t, y) = C 2x f 2y - l) if y > 1/2 
| Co* 1 G'.. ;>0 = (2.r — 1, 2y) if x > 1/2 and y < 1/2 
l Ojj 1 (x, y) — (2.r, 2yj otharwis 


By applying Algorithm (1) to this dynamical system 

with 

W = [(x, yd E M 2 : a < x < c, b < y < d} f where 

(c-, b) = (0, 0 J and (c, d'j = (l, lj and V 

= i(r, -yj G K z : r z E y z > 2} , / : (R z IR = , and A r is 

larger enough we get the fractal set S constructed by this 
algorithm defined by : 

S - {(i, y) e W : 6 R far all n < N). 

That is, the black point represents the set S. See figure (l.a) 
and we find the correlation dimension by this figure (l.b) 
show the plot of log C (if) against log(^) for the fractal. 
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correlation dimension is computed from the sloped linear 
portion of the plot determined from a least square fit. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


Figure (l.a) The Sierpinski triangle 

In the following log-log graph indicate the correlation 
dimension (D €or (S) = 1.5S6B) of the Sierpinski triangle 
generated by a computer Matlab program 


correlation dimention for fractal generated by Escape Time Algorithm 



Figure (l.b) 


w .n 

Example 2: Consider the IFS [IE£ where 

= (gxjjr).and u 2 <Jc,y) = (A + \^y) 

The attractor of the IFS is the fractal Cantor set C . The 
inverse transformations for the above system, which are as 
follows: 

Define / : M 5 , by 


r 

fb,yi =-j 



U.y) 
(x, yJ 


(3x, 3y) if x <\ 

(3x — 2,3y) t f x >\ 


In the same manner as in example (1), we calculate the 
correlation dimension of cantor set generated by Escape Time 
Algorithm, (jD ffi0r (C) = 0.6351)™ 


Fractal generated by Escape Time Algoritm 



Figure (2.a): The Cantor set 


correlation dimention for fractal generated by Escape Time Algorithm 
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Figure (2.b) 
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Consider 
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Example 3: 

[K 2 

lid-! Or, y) = ,w 2 Cx,y) = 

= (A 


IFS 

where 


4 2 - 4 

1 1 
4 ' 4 




!) 


The attractor of the IFS is the fractal T . The inverse 
transformations for the above system, which are as follows: 
Define / : K 2 , by 

f 

irjfH.iyy) = (4*, 4 y] 


1 


1 


fb,y) = 


fif 1 fc#y) = (4.v - 3, 
fi^Gr.y) = (2x-- t 2y-- 
^(x t y) - (4i,4y - 3j 


^j 5 L U, y) = (4r - 3,4_y - 3) 


7 *<5 d 

1 -r 
VI 

r-„. 
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In the same manner as in example (1), we calculate the 
correlation dimension of the fractal T generated by Escape 
Time Algorithm, (B ccr (jj = 1.4384). 




Fractal generated by Escape Time Algoritm 



Figure (3.a) 


correlation dimention for fractal generated by Escape Time Algorithm 



log 2 (s) 

Figure (3.b) 


Example 4: Consider the IFS 

A 

(R J , .. £rj.n fog , foj_, fog, fog, Ojj, fog 1 .. wlld'e 



The attractor of the IFS is the fractal Sierpinski Carpet Sc. 
The inverse transformations for the above system, which are 
as follows: 

Define / " M 2 K* 

by 


f(x,y) - 

r ttf 1 Cr,}0 = (3x—3y) 


c^ L (x f y) ~ (3* - 2,3y) 

if x > - , y < - 
1 a J a 

£^ 1 0r,y) = (3x. 3y - 2) 

if X <- ,y >- 

w7 l (x f y) = (3i - 2,3v - 2) 

( 

2 

if x>~,y >- 

fog L U,y) = (3-t, 3y - l) 

if x <- , - < v < - 

fog L Gr,y) = (3.r - 2, 3y - l) 

if x > - , - < y < 

j g 3 

fo 7 1 0r,y] = (3i - 1, 3y - 2] 

['/ - < X < - t y > 

J g g 

\ Wg 1 Of, yj = (3.r — 1, 3y) 



In the same manner as in example (1), we calculate the 
correlation dimension of sierpinski carpet generated by 
Escape Time Algorithm, (D €or (ScJ = 1.8924). 

Fractal generated by Escape Time Algoritm 

0. 

0. 

0. 

0. 

0. 

0 . 

0 . 

0 . 

0. 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x n 

Figure (4.a): The Sierpinski Carpet 


correlation dimention for fractal generated by Escape Time Algorithm 



Figure (4.b) 



Example 

^Cxvy) 

Or, y) 


5: Consider the IFS 

n a , 5 1 , 3 , a ) 

[-x —-y -h — ,-s F -y F — 1 

\2 B ^ 16 2 6^ 16/ 

fi 3 3 i 3 .in 

\2 B^ 16 2 B J 16/ 


where 


The attractor of the IFS is the fractal G. The inverse 
transformations for the above system, which are as follows: 
Define 
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/:i . 2 * US by fix.y) = 

= {* +■ y - Z’Y X + 3 ^ + i) if x < ’ 

(ys -1 Gf,y) = (x - y + ^ ,^x + - y) if x >~ 2 

In the same manner as in example (1), we calculate the 
correlation dimension of fractal G generated by Escape Time 
Algorithm. (D €or (.G) = 1.2990 )l 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 


Figure (5.a) 


correlation dimention for fractal generated by Escape Time Algorithm 



Figure (5.b) 

VI. Conclusions 

There are many approaches to compute the fractal dimension 
of an object. In this work we have presented an algorithm for 
computing correlation dimension D €cr of fractal generated 
by an Escape Time Algorithm (ETA) which is based on 
computing the correlation function. Computing the 
correlation function based on the selection of the Euclidean 
distance is presented and implemented the Matlab program 
listed in the Appendix for reducing the computational 
complexity of the Grassberger-procaccia Algorithm. 

APPENDIX 

Computer program for computing the correlation dimension 
of fractal generated by Escape Time Algorithm 
{-Graph the Shapes using ETA 

} 

clc; clear all; 

%a=l/3; b=0; c=0; d=l/3; 
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a=0; b=0; c= 1; d=l; 
numits = 4; r = 1; 

M=100; 

A1 =0; t=(); h=0; 
method = T; 
for p= 1 :M 
for q=l:M 

x = a + (c-a)*p/M; 
y = b + (d-b)*q/M; 
n = 0; 

while n <= numits 
n = n + 1; 

switch lower(method) 
case T 

%disp('Sierpinski Triangle') 
if (x>=0)&&(y>0.5) 
x=2*x; y=2*y-l; 
end; 

if (x>0.5 )&&(y<=0.5) 
x=2*x-l; y=2*y; 
end; 

if (x<=0.5) && (y<=0.5) 
x=2*x; y=2*y; 
end; 
case '2' 

%disp('Ex2') 
if (x<0.5) 

xl=x+y-l/2; yl=4/3*y-4/3*x+l/6; 
x=xl; y=yl; 
end; 

if (x>=0.5) 

xl=x-y+l/2; yl=4/3*x+4/3*y-7/6; 
x=xl; y=yl; 
end; 
case '3' 

%disp('Exll') 
if (x<=l/4) 
if (y<=l/4) 
x=4*x; y=4*y; 
end; 

if (y>=3/4) 
x=4*x; y=4*y-3; 
end; 

if (y>l/4)&&(y<3/4) 
x=2*x-l/2; y=2*y-l/2; 
end; 
else 

if (y<=l/4) 
x=4*x-3; y=4*y; 
end; 

if (y>=3/4) 
x=4*x-3; y=4*y-3; 
end; 

if (y>l/4)&&(y<3/4) 
x=2*x-l/2; y=2*y-l/2; 
end; 
end; 

case '4' 

%disp('Cantor Set') Ex 13 
if (x<=0.5) 
x=3*x; y=3*y; 
else 

x=3*x-2; y=3*y; 
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end; 
case '5' 

%disp('Exl6'); %disp('Sierpinski Carpet') 
if (x<=l/3) 
if (y<= 1/3) 
x=3*x; y=3*y; 
end; 

if (y>=2/3) 

x=3*x; y=3*y-2; 
end; 

if (y>l/3)&&(y<2/3) 
x=3*x; y=3*y-l; 
end; 
end; 

if x>=2/3 
if (y<= 1/3) 
x=3*x-2; y=3*y; 
end; 

if (y>=2/3) 
x=3*x-2; y=3*y-2; 
end; 

if (y>l/3)&&(y<2/3) 
x=3*x-2; y=3*y-l; 
end; 
end; 

if (x>= 1/3)&&(x<= 2/3) 
if (y<=2/3) 
x=3*x-l; y=3*y; 
end; 

if (y>=2/3) 
x=3*x-l; y=3*y-2; 
end; 

if (y>l/3)&&(y<2/3) 
x=3*x-2; y=3*y-l; 
end; 
end; 

otherwise 

disp( Unknown method. ) 
end 

if x A 2 + y A 2 > r 
t=t+l; 

P(t)=p; Q(t)=q; 

X(t)=x; Y(t)=y; 
n = numits+1; 

%fprintf('(%4.3f,%4.3f)=%4.3fVn',x,y,x+y); 

else 

h=h+l; 

Pl(h)=p; Ql(h)=q; 
end; 

if (x A 2 + y A 2) A 0.5 > R A 2 
A1 = A1 + 1; 
end; 
end; 
end; 
end; 

figure(l) 

plot(P(l :t),Q(l :t),'.blackVMarkerSize',5); 

xlabel('\itx_n') 

ylabel('\ity_n') 

title(' Fractal generated by Escape Time Algorithm '); 

{-Calculate correlation Dimensions 

} 

clear all;clc; 


Ntrans=1000; 

N_pts=3000; 
x0=0;y0=0; 
m=3; 

a=[0.5 0 0 0.5 0 0; 

0.5 0 0 0.5 0.5 0; 

0.5 0 0 0.5 0 0.5]; 
px=y0; py=x0; 
for j=l:Ntrans 

s=floor(m*rand+1); 
nx=a(s,l)*px+a(s,2)*py+a(s,5); 

ny=a(s,3)*px+a(s,4)*py+a(s,6); 

px=nx; py=ny; 
end; 
nx 
ny 

x=zeros(N_pts, 1); y=zeros(N_pts, 1); 
x(l)=nx; y(l)=ny; 
for j=l:N_pts-l 
s=floor(m*rand+1); 
x(j+l)=a(s,l)*x(j)+a(s,2)*y(j)+a(s,5); 
y0+l)-a(s,3)*x^)+a(s,4)*y(j)+a(s,6); 
end; 
x 

y 

figure(l); 
axis tight 

plot(x( 1: N_pts) ,y (1: N_pts),' .b' , 'MarkerSize' , 3); 

xlabel('\itx_n') 

ylabel('\ity_n') 

title(' Fractal generated by Escape Time Algorithm'); 
grid 

ED=sparse(N_pts,N_pts); 
k=(); 

for j=l: 100 %N_pts 
for i=j+1:100 %N_pts 

d=(x(i)-x(j)) A 2+(y(i)-y(j)) A 2; 

ED(i,j)=d; 

k=k+l; 

if mod(k,250000)==0 
fprintf('%d - %3.2f\n',k,d); 
end; 
end; 
end; 

ED=sqrt(ED); 

min_eps=double(min(min(ED+( 1000*ED==0)))); 
m_eps=double(max(max(ED))); 
max_eps=2 A ceil(log(m_eps)/log(2)); 
n_div=floor(double(log(max_eps/min_eps)/log(2))); 
n_eps=n_div+l; 

eps_vec=max_eps*2. A (-((l:n_eps)'-l)); 
Npairs=N_pts*(N_pts-1 )/2; 
c_eps=[]; 
for i=l:n_eps 
eps=eps_vec(i); 

N = (ED<eps & ED>0); 

S =double(sum(sum(N))); 
c_eps = [c_eps; S/Npairs]; 
end; 

omit_pts=3; 

kl=omit_pts+l; k2=n_eps-omit_pts; 
in_grid = kl:k2; 
xd=log(eps_vec)/log(2); 
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yd=log(c_eps)/log(2); 
xp=xd(in_grid); yp=yd(in_grid); 

[coeff,temp]=polyfit(xp,yp, 1); 

D_Cor=coeff(l); 
poly=D_Cor *xd+coeff(2); 

D_Cor 
figure(2); 
plot(xd,yd,'s- ! ); 
hold on 

plot(xd,poly,'r-'); 
axis tight 

plot([xd(kl ),xd(kl)], [-30,30] ; 

plot([xd(k2),xd(k2)], [-30,30] ,'k—') ; 
xlabel( log_2(\epsilon)'); 
ylabel('log_2(C(\epsilon))') ; 

title( {'correlation dimension for fractal generated by Escape 

Time Algorithm }); 

grid 
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